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I. INTRODUCTION 

As the characteristic properties of electronic circuits 
reaches quantum scales, new possibilities emerge for the 
realization of novel, quantum electronic devices which ex- 
ploit the quantum nature of solids. Here we show that 
by using quantum effects, the capacitance of electronic 
devices can be optimized to a great extent, which is a 
key t o the f urther miniaturization of electronic compo- 
nents.EIUlIll While in MOSFETs, for example, large gate 
capacitances are needed for operation at small gate volt- 
ages, their gate insulators are required to have an equiv- 
alent thickness of at least 1 nm to keep the gate currents 
at acceptably small values. To meet both requirements, 
high-K gate insulators, consisting for example of Hf02 
based compounds, have been introduced into the gate 
stack of MOSFETs used in microprocessors. 

As the sizes of the devices approach the nanometer 
length scale, quantum mechanical phemonena become 
increasingly relevant. The well-known classical models 
of semiconducting devices and the corresponding scaling 
laws break down. Ballistic transport processes, for exam- 
ple, gain importance, as does the statistical distribution 
of the dopant atoms in the small drain-source channels. 

Quantum effects also alter the capacitance of large, 
macroscopic capacitors that reach the quantum regime 
because their dielectrics become thin. As we show, for 
such devices the capacitance comprises not only the well- 
known geometric capacitance, but four more capacitance 
terms. These originate from the kinetic energy of the 
electrons, from their exchange energy, their correlation 
energy, and the electron-phonon coupling energy. 

In this paper, we present the general equation for the 
capacitance of capacitors that considers the quantum me- 
chanical energies of the electronic systems involved. This 
equation applies to two-plate capacitors, also in case their 
dielectric is ultrathin. It is shown that for given and 
A the capacitance of capacitors reaching the quantum 
regime may differ considerably from the well-known value 



In fact, the capacitance may be even negative. Here cq is 
the dielectric constant of the vacuum, is the dielectric 
constant of the dielectric material between the two elec- 
trodes of area A, and d is the thickness of the dielectric. 

It is specifically shown that capacitors can be fabri- 
cated that do not follow Eq. ([ij. By altering the con- 
ducting materials of the electrodes and by varying the 
geometry of the capacitor in novel ways the capacitance 
of such capacitors can be greatly enhanced or lowered as 
compared to Eq. ([T|. Electrode materials of interest in- 
clude materials with strong electronic correlations such as 
transition metal oxides, two-dimensional electron gases 
as generated at interfaces in oxides or semiconductors, 
metals with small carrier densities, organic conductors, 
and materials such as graphene. We point out that these 
quantum phenomena and their technical implementation 
provide a possible strategy to control large or small ca- 
pacitances by the electronic properties of the conducting 
electrodes — as opposed to the standard search for low 
or high-K dielectric^^. 

It is common knowledge that classical electrostatic^^ 
determines the capacitances. In the textbook case of a 
capacitor with two electrodes, the capacitance is charac- 
terized by 

Q = CV (2) 

where Q is the charge on a capacitor electrode, C is 
the capacitance, and V is the voltage between the elec- 
trodes. In more general terms, for multiple conductors 
{i — 1,2, . . . ,n), the capacitance is characterized by the 
corresponding static response relations Qi — J^j CijVj, 
where Vt and Qi are the potential and the total charge 
of conductor i, respectively. The coefficients de- 
fine the capacitance matrix. According to the classi- 
cal approach, the coefficients of a capacitor consisting 
of two infinitely extended parallel metallic plates are 
Cii = C22 = -C12 = -C21 = C (Eq. (0). Here, fringe 
field effects have been neglected. In this classical descrip- 
tion, the capacitance is not infiuenced by the electronic 
properties of the conducting plates. For a two-plate ca- 
pacitor, the geometric capacitance is entirely determined 
by the effective distance d* = d/er- 
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It is apparent that Eq. ([T| is not very useful for the 
design of ultra-small devices, where, for example, single 
electron effects play an essential role (see, e.g., Ref. [T6j. 
The calculation of the electronic properties of such de- 
vices is typically based on the direct calculation of quan- 
tum mechanical energies. 

Quantum mechanical effects, however, are also of im- 
portance for large devices with macroscopic lateral di- 
mensions. The conductors of electrodes can screen, for 
example, the electric field only over non-vanishing lengths 
such as the Thomas Fermi screening length. This effect 
obviously enhances the effective dielectric thickness d*. 
The enhancement of d* has first been observed by MeacP, 
confirmed by Hebard et alW for a known thickness of the 
dielectric, and deducecP^ from the charge distribution 
profile in the electrodes which is controlled by the kinetic 
energy of the charge carriers (see the detailed discussion 
and evaluations in Refs.[9j. A framework to calculate the 
capacitance coefficients Cij that considers the field pen- 
etration into the conductor has been formulated by Biit- 
tiker.Ii^He evaluated the space-dependent Lindhard func- 
tion. The field penetration into the capacitor electrodes 
introduces corrections to the geometrical capacitanceP^ 
in the case of space independent density of states (DOS) 
right up to the surface A/C = {d/e^ + A]^ + ^2 )/eo 
where \\ 2 are the Thomas-Fermi screening lengths 
of the two metal plates. For a capacitor with two- 
dimensional metallic plates, A^"'"^'' -I- A2"''^^ is replaceci^^l^ 
by 03/2 where as = 47reo/i^/(me^) is the Bohr radius. 

Lang and Kohi^Sl introduced the center of mass posi- 
tions a^g ' ' of the induced surface charge density at the 
electrodes and related the distance between the center of 
mass positions to the capacitance of a parallel plate ca- 
pacitor in the vacuum (e,. = 1): A/C = {x^q^ — x'q^)/^^). 

(i) 

They calculated Xq within density-functional theory 
(DFT^ for a metal with a uniform distribution of pos- 
itive background charge (jellium model). They deter- 
mined the center of mass positions to be outside of the 

surface of the electrode (given by the edge of the pos- 

(1 2) 

itive background charge at xl ) so that the capaci- 
tance is increased with respect to its geometric value 
60^/(2:^^'' — x^^^). The information on the electrodes is 

encoded in the distances fe^^'^-* — Xq^'^'' — x^^'^-* which 
generate interface capacitances serially connected to the 
geometric capacitance. 

For a macroscopic d*, the corrections arising from the 
Thomas-Fermi lengths or are tiny, because they 

represent small, quantum mechanical length scales. How- 
ever, exceeding these effects described, quantum phenom- 
ena are of fundamental relevance in controlling the capac- 
itances of electronic circuits as soon as the effective thick- 
ness of their dielectric layers d* approaches quantum me- 
chanical lengths such as erUB- Current devices like mod- 
ern MOSFETb are close to this transition, and quantum 
phenomena need to be considered in their design. As 
we show, this even holds for capacitors of macroscopic 



lateral dimensions and macroscopic numbers of charge 
carriers, if their dielectric layers are sufficiently thin. It 
is hereby particularly interesting that also correlation ef- 
fects in the electron systems of the conducting electrodes 
may alter the capacitances in unique ways. We therefore 
propose to use these quantum phenomena as a possible 
tool to optimize devices such as MOSFETs. Desired large 
or small capacitances can potentially be attained by tun- 
ing the electronic properties of the conducting electrodes 
— in addition to the standard approach of using low- or 
high-K dielectrics and small or large barrier thicknesses. 

II. THE COMPOSITION OF THE 
CAPACITANCE OF MACROSCOPIC 
CAPACITORS 

To determine the equation that describes the capaci- 
tance of capacitors, we start by considering the energy 
functional for the electronic ground state of a given ca- 
pacitor. For clarity, we restrict our evaluation to capaci- 
tors with two parallel plates and neglect effects arising 
from finite temperatures and from fringe fields. The 
plates are considered to be so far apart that their elec- 
tronic wave functions do not overlap. Tunneling and 
exchange contributions between the plates are therefore 
negligible.!^ The ground state energy as a functional of 
the electronic densities may be decomposed asJ^ 

E = E'h + ^ -Ekin/i + ^ Ex,i + ^ Ecs + ^ Ecxt.i (3) 

i i i i 

The sums are used to describe the energies of the two 
possibly different metallic plates. In capacitors for which 
surfaces are relevant, the electronic density nj(r) is a 
function of the space coordinate r. The energy E[ni{r)] 
is then a functional of the densities on the two electrodes. 

The Hartree term En in Eq. ([3} includes all direct 
Coulomb interactions between the electrons of the con- 
duction bands; these long-range Coulomb potentials also 
include the electrostatic interaction between the plates. 
The other energies, which are of pure electronic origin, 
are not caused by inter-plate interactions, but rather are 
local attributes of each plate. These energies are the ki- 
netic energy of the electrons i?kin.i, the exchange or Fock 
energy Ey^ i, and the correlation energy E^^i, which rep- 
resents the electron-electron interaction energies beyond 
the Hartree-Fock description. 

The external energy i?cxt,i is the energy of the conduc- 
tion electrons in an external potential; it most promi- 
nently includes the interaction energy between the elec- 
trons and the background charge that is caused by the 
nuclei. Because this Coulomb interaction between elec- 
trons and nuclei is long-ranged, it diverges with increas- 
ing area A; this divergence is canceled by a similar diverg- 
ing term of opposite sign, which results from the electron- 
electron interaction in the Hartree contribution.^S Inho- 
mogeneities such as surface potentials or potential wells 
are also introduced through -Ecxt.i- 
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To treat the electronic exchange effects, all electron- 
electron interactions, including those with the electrons 
of the ionic core states, need to be contained in the 
Hartree-Fock and correlation terms. In effective models, 
however, the exchange contribution of ionic core states 
with conduction electrons is often neglected, an approx- 
imation that allows to include the occupation of these 
electronic core states in a pseudopotential or in a suit- 
ably defined ionic core charge. In those cases, it is useful 
to introduce an effective dielectric constant €es,i into the 
electronic interaction terms which accounts for the po- 
larizability of the ionic cores. 

In case the electron systems on the two parallel plates 
are homogeneous, i.e. if they don't reside on lattices,'22l 
En + Ecxt,i — Q^d/{2eoerA). In such a situation, d is 
given by the distance between the boundaries of the uni- 
form positive background charge of the electrodes. This 
energy term matches the standard equation for the en- 
ergy of the capacitor. 




FIG. 1: Equivalent circuit diagram of a two-plate capacitor 
according to Eq. ([gJ. The capacitor consists of a serial con- 
nection of the geometric capacitance and the quantum capac- 
itances of each of the two plates. The capacitance results 
from the Hartree energy (oc 1/Cgeom), the kinetic energy 
(oc 1/Ckin), the exchange energy (oc 1/Cx), the correlation 
energy (oc 1/Cc), and the short range interaction energy of 
charge carriers with all "non charge-carrier" (ncc) degrees of 
freedom (oc 1/Cei-ncc)- Positive and negative capacitances 
are drawn in both blue and red, respectively. The capaci- 
tance Coi-ncc is presented in blue and red, because the sign 
of the capacitance depends on the microscopic model for the 
coupling of the ncc degrees freedom to the charge carriers; for 
phonons it is expected to be negative. 



^'Coul — Q^/(2Cgcom) (4) 

with the geometric capacitance Cgcom — ^o^rA/d. Yet, 
also the other three terms in the energy functional gener- 
ate capacitances. These capacitances result from distinct 
quantum phenomena. 

To derive the total capacitance C for Q — > 0, we ex- 
pand the energy functional to second order in Q. If the 
two plates differ in their physical properties, also a lin- 
ear contribution in Q arises, which originates from the 
difference in the work functions of the electrode materi- 
als. The work function defined as the minimum energy 
necessary to extract an electron, is composed of the chem- 
ical potential of the infinite system (relative to the elec- 
trostatic potential in this system) and the electrostatic 
potential energy difference across the respective surface 
of the electrode (for details and the evaluation see, e.g., 
Refs. El] and [T5]l . The contribution to the work func- 
tion from the chemical potential is generated by all en- 
ergy terms. The electrostatic potential difference across 
the surface, the so-called electrostatic dipole barrier, is 
caused by the external energy. A difference A<I> of the 
work functions of the surfaces of the opposing electrodes 
generates a contact voltage, so that associated surface 
charges ±Qs are induced. This electronic charge recon- 
struction compensates the work function mismatch. The 
differential capacitance can be evaluated by expanding 
the energy functional around the charge-reconstructed 
states of the electrodes. 

Focussing first on the capacitance 

C = Q/V (5) 

in linear response^ i.e., for V ^ Q (or, for finite 

= -CA$/e, 01J22] c = (Q- Qs)/V), we collect the 
quadratic terms and identify the inverse capacitances ac- 



cording to the coefficients of /2: 

1/C=l/Cgoom + 1/ Ckin.i + l/C'x.z 

i i 

+ ^l/Cc,.-H^l/Col_„cc,, (6) 

i i 

Eq. ([ij is replaced by this general equation, which reveals 
that the total capacitance is caused by a series connection 
of the capacitances that result from the additive terms in 
the energy functional (Fig. [iJ. Although Eq. (jlj seems to 
imply that the total capacitance is always smaller than 
the geometrical capacitance, the total capacitance can, 
in fact, well exceed the geometrical one, because Cx and 
Cc may have large negative values. 

Note that the external energy i?oxt,i consists of a long 
range Coulomb term, which compensates the diverging 
term of the purely electronic Hartree energy, and a re- 
maining short range contribution. The inverse geomet- 
ric capacitance 1 / Cgoom is generically defined as the Q^- 
coefRcient from the Hartree energy and the static long- 
range part of the external energy. The short range part 
of i?oxt,i contribution includes the interaction between 
the charge carriers and all "non-charge carrier degrees of 
freedom" (ncc) such as phonons or local potentials. Ex- 
amples for the latter are the static potentials of the ionic 
cores and the confining potentials in semiconductor het- 
erojunctions. These contributions generate the l/Coi_ncc 
term in Eq. ([6}. In this paper we will not investigate 
Cci-ncc in more detail. 

Eventually, to obtain the voltage dependent capaci- 
tance C{V) beyond the linear response limit of Eqs. ([5} 
and ([6} , one takes the Legendre transform E of the energy 
E{Q) with respect to Q and identifies the capacitance 
C{V) = Q{V)/V from E{V) = E{Q{V)) - QV. 

For the sake of clarity, we explicitly note that the value 
of the differential capacitance for specific values of V, Q, 
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or the density of charge carriers, n, is given by the ratio 
dQ/dV, where dQ is the amount of charge that a volt- 
age change dV adds to a plate of the capacitor, which is 
characterized by V, Q, or n. We also note that the volt- 
age V between the two plates, i.e., the difference between 
the electrochemical potentials of the plates, does in gen- 
eral not equal the difference of the electrical potentials 
of the electrodes, which defines the electrical voltage Vc- 
Whereas the voltage V is the voltage that is measured 
by standard voltmeters and, for example, is given by the 
output of standard voltage sources, the electrical volt- 
age Vc is the electrical potential difference between the 
plates, the gradient of which is the electric field between 
the plates. The electrochemical potential difference con- 
sists of the electrical voltage and the difference of the 
chemical potentials of the electrodes. In standard capac- 
itors, the shifts of the chemical potential upon charging 
are negligible. However, if the chemical potential differ- 
ence becomes significant, as is the case when the charging 
of the capacitors cause a non-negligible shift of the chem- 
ical potentials, the measured voltage will differ from the 
electric potential difference. 



III. CAPACITORS WITH ELECTRODES 
COMPRISING ELECTRON SYSTEMS WITH 
UNIFORM BACKGROUND 

Electron gases with uniform, positively charged back- 
ground (jellium model) serve as model systems which al- 
low to evaluate thermodynamic electronic properties in 
2DE2] and 3DES1 without the complications arising from 
the coupling to lattice degrees of freedom. In the fol- 
lowing we will calculate the capacitance for 2D or 3D 
electron gases in the jellium model to make several of 
our statements on the quantum mechanical properties of 
capacitors more lucid. The ground state of the electron 
gases in a configuration with two unbiased but electro- 
statically coupled sheets, each comprising a dilute elec- 
tron gas, has been investigated intensely. In particular, 
it has been laid open that the exchange term ^x,i can 
drive the double-layer el ectron g as into ferromagnetic and 
charge separated states.'22EIIlir 

A spontaneous spin polarization of a dilute electron gas 
was already suggested by Felix Bloch!^ and has been in- 
cluded in several recent investigations of double layer sys- 
tems. An exchange driven instability towards a state in 
which all electrons occupy a single layer, i.e., a state with 
a maximally charged capacitor has been proposecP2] but 
it has been shown that the ground state in an unbiased 
double layer system with E = Ecou\ + J2i ^kin,i + X]i 
is of a different natureP3] it is a coherent quantum me- 
chanical superposition of the two states that each carry 
the electrons on either of the two layers. This ground 
state has been introduced and predicted for bilayer quan- 
tum Hall systems. 12^121] n jg characterized by a spon- 
taneous interlayer phase coherence (SILC) of two lay- 
ers which are each entirely isolated except for their in- 



terlayer Coulomb interaction. The mean value of the 
charge in the SILC state is zero on each of the two lay- 
ers. The Coulomb energy i?coui therefore vanishes in 
this state. The possible magnetic, charge separated, and 
SILC states have been categorized in Hartree-Fock theory 
at zero magnetic field, with and without applied bias.ISS 
A finite bias may support charge separated states. 

Spin polarized and SILC states offer exciting possibil- 
ities to explore and use quantum capacitances. Never- 
theless, in the rest of this article we will focus on un- 
polarized states without interlayer coherence. While the 
SILC states will be the true ground states, coherency 
will be lost in capacities operating with finite electric 
fields in noisy electronic circuits at room temperature. 
The states with charged electrodes may be generated 
through an applied voltage in cases where the differential 
capacitance C for Q ^ is negative (see the following 
subsections). The states with charged electrodes are ex- 
pected to persist up to high temperatures of the order 
of |£'coui + J2i -Skin.i + J2t E^.t + J2t Ec,t\/kB. In com- 
pounds with strongly electronic correlations, the correla- 
tions may favor either the SILC states or the states with 
charged electrodes. We consider the unpolarized homo- 
geneous electron gas on two nanoscopically close plates 
as a "model capacitor" and refer to possibly relevant, cor- 
related electron systems in the following chapters. 

We conceive that practical capacitors with negative, 
total capacitance can eventually be realized, in particu- 
lar by using strongly correlated systems. The existence of 
devices with negative total capacitance has already been 
addressed in the literature, concerning electronic circuits 
(see, e.g, Ref. |39| . electron injection through interfaces 
(see, e.g, Ref. HOj, and structures involving ferroelectrics 
(see, e.g, Refs. HTl H2]) . Here, we consider a far more con- 
ventional system, standard capacitors, such as plate ca- 
pacitors comprising conventional dielectric materials and 
operated in steady state. 



A. Electrodes with Two-Dimensional, 
Homogeneous Electron Systems 

The capacitance C for a capacitor with two- 
dimensional electron systems on each plate is determined 
by the inverse capacitances given in Eq. To shed 

light on the physical origin of the rhs of Eq. ([6| and to 
present a method to approximately calculate these terms, 
we introduce the thermodynamic compressibility x: of an 
electron system. As follows from Seitz's theorenP^, 
is given by the first derivative of the chemical potential 
/i or by the second derivative of the total energy E, both 
taken with respect to the electronic density n — N/A: 
— n^dn/dn — n? d'^ {E / A) / dn^ , which is a standard 
thermodynamic relation in the limit of zero temperature. 
For 2D-systems we rewrite Eq. (|6| in terms of the inverse 
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compressibility of a 2D conductorP 



(7) 



The inverse kinetic and exchange capacitances can be 
determined from Eq. ([7|, using the well-known expres- 
sions for the energy functional of a 2D homogeneous elec- 
tron system (see, e.g., Ref . [Mil : 
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47reov4/a 



m* /m 
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2 1 1 1 



(8) 
(9) 
(10) 



where qb — '^ire^h^ / [me^) ~ O.SA is the bare Bohr ra- 
dius, m* is the effective electronic mass and eoff,i the ef- 
fective dielectric constant in plate i. This constant char- 
acterizes the screening in the metallic 2D plate that does 
not arise from the screening by the conduction band elec- 
trons themselves. The band electron "self-screening" is 
included in the correlation term Cc discussed below. It is 
notoriously difficult to estimate eofi,i. For very dilute 2D 
electron systems in GaAs/AlGaAs heterostructures that 
may have nominal densities of 10^^ cm~^, Coff is given 
approximately by the dielectric constant of the AlGaAs 
layer.Ssl The effective band mass m* can be more eas- 
ily determined. Conveniently it is chosen such that the 
resulting DOS matches the DOS pf'^\£F) of the con- 
ductor. A summation over bands is used in the DOS 
in case several conduction bands are present or a band 
degeneracy exists.l^ Such a band (degeneracy) factor de- 
creases the inverse kinetic capacitance and, consequently, 
increases the total capacitance C. 

The kinetic capacitance C^™], which has been intro- 
duced by LuryiPS as "quantum capacitance", is invariably 
positive. The exchange capacitance A/Cx, however, is 
always negative, irrespective of the dimensionality of the 
electron gas. The exchange capacitance varies inversely 
with the square root of the electronic density. Corre- 
spondingly, based on Eq. ([6| we recognize that capacitors 
with negative capacitances can be built. The negative 
capacitance can be obtained, for example, by fabricating 
capacitors with electrode materials that have a very low 
density of states. A negative C implies that the plates 
of the capacitor self-charge once they are connected. In 
these systems, the inverse exchange A/Cy^ plus correla- 
tion A/Cc capacitances dominate the inverse kinetic ca- 
pacitances A/Ckin plus the geometrical capacitance, so 
that the total capacitance C becomes negative. 

In a homogeneous electron gas, electron-phonon contri- 
butions are not present by definition. Correspondingly, 
A/Cci-ncc in Eq. ^ is disregarded here. 

The last term to be considered in A/C of Eq. ([6} is the 
correlation term. If present, Cc is particularly important. 



and obviously even more so in those cases in which A/C^ 
and A/Ckin almost cancel each other. 

To estimate the correlation capacitance of the 2D ho- 
mogeneous electron gas, the parameter is introduced. 
It is a dimensionless number that characterizes the inter- 
particle distance: 



1 
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(11) 



Tanatar and Ceperlej^^U have cast the r^-dependence of 
the correlation energy into an approximate functional 
form which fits their Monte Carlo data. With the en- 
ergy functional of Ref. [Ml the inverse capacity is identi- 
fied by using the correlation contribution to Eq. Q. The 
exchange and correlation capacitances are found to equal 



and 
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ao / 3 3,2: , ,a; z.\ 

- -T^"+(l + T«)- + (l + «)(^--)K 
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respectively. Here, the polynomial functions u, x, y, z are 
defined as: 



y[rs] = 1 + airi + aar^ -|- a^ri 

r 1 1 * 1 3 
x[rs\ = -airi + -a2rs + vase's 



4 
5 

32 



airi 



-a2rs 



21 
32 



flag's 



(14a) 
(14b) 
(14c) 

(14d) 



and the coefficients a; ard^ 



ao = -0.3568, oi = 1.1300, 03 = 0.9052, 03 = 0.4165 

The inverse capacity of each of the plates contributes a 
term A/C^l''^ to the total inverse capacity. For large r^, 

i.e., for dilute, homogeneous electron systems, A/c'^^^ oc 
— Ts cx —\j^n. This implies that in systems with very 
few carriers, C is driven by Cx and Cc to small, negative 
values (Fig. [2]). 

Eqs. ^ and (|9|, ([iT]) and ([l3| reveal that the total 
capacitance can differ considerably from the geometri- 
cal capacitance. It is only for the case that Cgcom/^ 
is much smaller than |>fj| (en^)^ that the total capaci- 
tance is bound to be given by the geometrical capaci- 
tance. The strong dependence of C on is illustrated 
in Fig. [2] for a symmetric capacitor. As expected from 
the previous discussion, with increasing the capaci- 
tance diverges and becomes negative when the negative 
terms 2(A/Cx -I- AjCo) compensate the positive terms 
(2A/Ckin + A/Cgoom) at Tg — tq. The value at which 
the capacitance changes sign is higher than the one at 
which the compressibility becomes negative. We assign to 
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FIG. 2: Dependence of the capacitance of an ideal two- 
dimensional capacitor (a parallel plate capacitor with two 2D- 
electrodes) on the carrier spacing. The capacitance is a func- 
tion of the electron density on the electrodes, parameterized 

by the dimensionless interparticle distance = (l/agTrn)^ 
in bare atomic units. In the calculation presented, the effec- 
tive mass is taken to be m*/m = 1, the effective dielectric 
constant in the electrodes to be ecff,i ~ 1, and the effective 
distance to be d* = as. The dotted line shows the capaci- 
tance of the conventional, classical capacitor Co = Cgoom. 



the latter the dimensionless length scale r^, and rg > Tc- 
In fact, the transition into a state of negative electronic 
compressibility has already been experimentally observed 
at interface e lectron gases in Si-MOSFETs and III-V het- 
erostructured^^E^IIZl_ The observation of a negative k 
does not necessarily contradict the thermodynamic sta- 
bility criterium of a positive total compressibility. The 
latter condition applies to the sum of the compressibil- 
ities of all electronic and ionic subsystems. Even in 
thermodynamic equilibrium, the compressibility x of the 
electronic subsystem on one or on both of the electrodes 
may be negativ^^U if it is compensated by other contri- 
butions to the total compressibility. 

It has to be expected that the electronic states on 
the plates change qualitatively at the transition into 
the regime of negative electronic compressibilitjH^] at Tc. 
The new ground state is characterized by exponentially 
damped charge density modulations (cf. [l9]). A further 
transition into a self-charged state takes place at rg where 
the capacity becomes negative. 

As shown by Fig. [2] the capacitance exceeds the clas- 
sical capacitance (C > Cgoom) for Tc < < rg. This 
result relies on the assumption that, also at these car- 
rier densities, the electron system is homogeneous. For 
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m*/m ^ 

FIG. 3: Contour plot, showing for an ideal two-dimensional 
symmetric capacitor the lines of constant capacitance. The 
capacitance C^^ is a function of the effective mass m* and 
of the effective dielectric constant eoff of the electrode ma- 
terials (assumed to be equal on both electrodes). Here, the 

dimensionless interparticle distance rs = (l/a%TTn)^ in bare 
atomic units is taken to be rs = 7 and the effective distance 
of the capacitor plates to be d* = d/er — 1.23aB, where is 
the dielectric constant of the dielectric material between the 
plates. The white line traces the values of m*/m versus ecs 
at which the capacitance diverges (cf. Fig. |2|. 



non-homogeneous electron systems, the phase transition 
into an electronic state with negative compressibility has 
not been investigated rigorously. It cannot be excluded 
that for Ts > Tc the capacitance behaves differently than 
shown by Fig. |2] In particular, a metal-insulator transi- 
tioii^with a spatially inhomogeneous insulating phasd^ 
has been observed which is probably controlled by strong 
disorder. Also, a phase transition into a charge density 
wave state seems possible.!^ The functional properties 
of the capacitance nevertheless reflect the energy depen- 
dence of the electronic state on the charge carrier density 
of the electrodes. 

For the ideal capacitor with two-dimensional elec- 
trodes, the kinetic term has a negligible influence on the 
capacitance for large density of states or, equivalently 
for large effective mass (see Fig. [sj. This applies for 
P^^^H^-f) > 2eoer/{e'^d), i.e., for m*/m > eras I {2d). 

Leaving the model capacitor, the description of a real- 
istic quasi-2D electronic system has to consider also the 
finite thickness of the quantum well which hosts the 2D 
electron system. As discussed, the Hartree term con- 
tains such corrections which are included in d* . Also 
the exchange and correlation terms depend on the finite 
thickness which may be considered by a "form factor" 
that acts as a prefactor for the respective energies. Such 
a form factor has been introduced by SterrP for the ex- 



7 



change term and has been applied by Eisenstein et aZJ^ 
in the evaluation of their experiments. This finite thick- 
ness correction does not alter in a qualitative way the 
dependence of the capacitances as given by Eq. [T2| 

Furthermore, also the infiuence of the charge carrier 
density on the band energy (band bending) has not been 
considered in the model capacitor model. This effect de- 
pends on the details of the electronic structure of the 
electrode material. A detailed analysis is given in Ref.fiSl 



B. Electrodes with Three-Dimensional, Metallic 
Electron Systems 

Quantum corrections to the capacitance that result 
from the kinetic energy of the additional carriers of a 
charged 3D capacitor, have been investigated in Refs. jO] 
I12I55I If the corresponding states are confined to the sur- 
face, i.e., if they form subbands of surface bound states, 
they are two-dimensional states. Their existence requires 
that their energies are well separated on the temperature 
scale ksT and that the inelastic scattering between states 
of these subbands is negligible. In that case the evalua- 
tion of Sec. [Ill A] applies. In this subsection we consider 
the case that the electronic states are not bound to the 
surface of the electrode. 

Also in 3D systems, the geometric capacitance is given 
by Eq. (§, 

47reoA/Cgeom = 47rd* (15) 

except that the charge does not completely reside in the 
surface planes of the electrodes. A simple approach for 
3D systems assumes the charge density to decay into the 
bulk with a length scale of the order of the screening 
length A. The effective distance between the capacitor 
plates is thereby slightly enhanced by 2A, because the 
static Coulomb energy is determined by the charge dis- 
tribution which is infiuenced by the bound ary condition 
of the electron wave function at the surface.l^SllSnilMlII] 
The results of such an approach have to be contrasted 
with those of Lang and Kohn.lSl They include the for- 
mation of a surface dipole barrier that is caused by 
the "spilling out" of the mobile electrons into the vac- 
uum. Already in the uncharged state this dipole layer 
forms a non-homogenous charge profile. This charge 
profile modifies the distribution of the induced screen- 
ing charge for finite voltage so that the effective distance 
between the electrodes is reduced rather than enhanced 
by 2A, assuming that the difference between the edge 
of the uniform background charges defines the distance 
d*. Whereas the width of the screening charge density 
scales approximately with the Thomas-Fermi screening 
legnthP^^ — apart from the characteristic Friedel oscilla- 
tions in the interior of the metal — the center of mass 
position of the induced charge is not easy to determine. 
Its calculation requires a fully self-consistent treatment 
that includes the exchange and correlation terms. 



To elucidate the behavior of capacitors with metallic 
electrodes that contain 3D electron gases,'2Sl we consider 
in the following a model capacitor. In this capacitor, 
the 3D bulk compressibility of the electronic systems is 
taken to be constant up to the electrode surfaces; con- 
sequently the formation of a surface dipole layer is not 
included. With this approximation, the precise distance 
between the mirror planes of the electrodes cannot be 
identified, as in Ref. [T3]for the jellium-model electrodes 
separated by vacuum. However a comparison of the qual- 
itative dependence of the effective distance on exchange 
and correlation effects is of considerable interest and will 
be presented below. This allows to introduce an effective 
model which captures an important part of the physics 
and is much easier to treat than the fully self-consistent 
evaluation of Ref. [l3] It may seem that the approach, 
which is presented below and builds on the evaluation 
of the density-density correlation, cannot account for the 
displacement of the center of mass of the induced charge 
density as introduced in Ref. [T31 The density-density 
correlation in the uncharged state of the electrodes is, 
however, related via the fiuctuation-dissipation theorem 
to the displacement of the induced charge for finite, pos- 
itive compressibility. 

In three dimensions, the parameter is defined by 



and presents, like in the 2D systems, the dimension- 
less length scale of the carrier spacing. We assume that 
the electrodes in the uncharged state consist of a homo- 
geneous electron system with a uniform positive back- 
ground which compensates the electron charge (jellium 
model). Correspondingly, and the energy E are not 
functional of ni[x] but just functions of the average elec- 
tronic densities as in Sec. |III A) In a microscopic evalu- 
ation of the work function of the electrodes, this approx- 
imation would miss the electrostatic potential across the 
metal surface, because the homogeneous electron system 
(up to the edge of the positive background) does not in- 
clude the surface dipole barrier. However, these terms 
should be included phenomenologically in the difference 
of the work functions for the two electrodes in the capac- 
itance. For electrode plates of the same material and the 
same crystallographic orientation of their inside parallel 
faces, these terms cancel. 

With the 3D electronic density n = N/V, the 
compressibility relation is: >c~^ = n^dfi/dn = 
m? (P{E/V)/d'n?. For a capacitor with 3D plates, the 
screening of the charge causes the charge density to be 
non-uniform perpendicular to the plates on the scale of 
the screening lengths A^. With the above introduced 
assumption of constant compressibility, the relation be- 
tween the capacitances and compressibilities has the ap- 
proximate form 

= E 4rl^ - E (17) 

gcom ^ ^2\^ Qj^. Ai(eni)2 
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In case the density of states (DOS) p^'^'^\eF) (the DOS 
at the Fermi energy ep for both spin directions) is spa- 
tially uniform and if the screening follows the Thomas- 

('TP) — ^ 

Fermi screening of a free electron gas with A, = 
tQ^e^ pf^\eF), Eq. ( [l7| l is equivalent to the relation 
A/C7(3D) _ ACg-i^ = Y.^ Af^^ given by BiittikeiCll. 
The relevant surface volume of the electrode i is given by 
the surface area times the screening length A^. 

Eq. (171 yields the value of the inverse capacitance. 



The kinetic term in the energy functional generates for 
each electrode a capacitance, the inverse value of which 



as n. 



47ren 



\i m*/m 



(18) 

A straightforward conversion to the Thomas-Fermi ap- 
proach, which also accounts for the polarizability of the 
underlying lattice and of the ionic cores through a cor- 
responding effective dielectric constant e^g^i (see, e.g., 
Ref. [13}, identifies: 



A,; 



(TF) 



(19) 



where the bulk Thomas Fermi screening length is 



,(TF) 



1 




(20) 

Eq. ( 18 1 Jointly with Eqs. ( 20 1 and (191, yields the simple 
relationlSl 



47rA, 



(TF) 



(21) 




FIG. 4: Dependence of the capacitance of a three- 
dimensional, parallel plate capacitor on the carrier spac- 
ing, according to the serial connection of the capacitances 
of Eqs. (pSl), (irsl), (|23l and (ISil). The capacitance is a 



function of the electron density on the surface of the elec- 
trodes, parameterized by the dimensionless interparticle dis- 
tance rs = (3/47rnaa)5. In the calculation shown, the effec- 
tive mass is m* jm = 1, the effective dielectric constant in the 
plates is tc9L,i = 1, and their effective distance is d* — as- 
The dotted line displays the capacitance of the correspond- 



ing classical capacitor Co 
to a capacitor with ecff,i - 
capacitor with eoff,i: ~ 30. 



The dashed line refers 



2, and the dot-dashed line to a 



Here we have used the density of states (DOS) relation 
for an electron gas: 



(^^) = (zi) 



3 - i Anenns m* 



1 . 3 47reo m 



TT 



m 
(22) 

We derived the relation Eq. (21 1 from the (constant) com- 



pressibility of the kinetic energy term, whereas the stan- 
dard derivation determines first the induced charge den- 
sity in Thomas-Fermi approximation. Both schemes are 
obviously equivalent if the screening length is identified 
as A^). 

The quantum kinetic term (Eq. ([l8|) depends on the 
DOS and always lowers the capacitance as compared to 



its classical value. As shown by Eqs. (|T8| and (20 1, mul 



tiple or degenerate conduction bands in the DOS increase 
(7(3D)_ Conversely, small effective masses m* of the car- 
riers reduce C^^^^ (cf. Fig. [HJl. 

Coulomb interactions beyond the classical approxi- 
mation qualitatively alter the properties of the capaci- 
tance, because exchange and correlation terms introduce 
in Eq. (IgIi negative inverse capacitances. With increasing 



electronic correlations the capacitance therefore grows to 
jump at Ts > tq to negative values (see Fig. [ij. 

The exchange contribution to the capacitance results 
from the standard exchange energy functionaff^ (using 

Eq. my. 



47reoA/a 



(3D) ^ 



1 1 

(Q'^)^ A,eeff,.n| 
/167r\5 n a'j^ 



AiEeft,; 



(23) 



The correlation contribution is derived from the ap- 
proximate functional form of the 3D energy fit to Monte 
Carlo calculations as given by CeperleypS 



47ra|r3 



3Aiecff,i f[rs 



,9Fs 



9FS 



f[r 



(24) 
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FIG. 5: Dependence of the capacitance of a three- 
dimensional, parallel-plate capacitor as a function of carrier 
spacing. The effective distance between the plates used in the 
calculation is d* = as and m*/m = 5. The effective dielectric 
constants are equal for both electrodes with eeff.i = 1,2, and 
30 for the solid, dashed, and dot-dashed lines, respectively. 
The dotted line shows the capacitance of the corresponding 
classical capacitor Co = Cgeom- 



FIG. 6: Contour lines of equal capacitance for the three- 
dimensional model capacitor. The capacitance is a func- 
tion of the effective mass m* and of the effective dielectric 
constant e^f! of the electrode materials. Here, the dimension- 
less interparticle distance rs — (3/ATvna%)^ in bare atomic 
units is taken to be = 7 and the effective distance of the 
capacitor plates to be d* = d/er = 1.23 as. The gray line 
traces the values of m* /ra versus e^s at which the capaci- 
tance diverges. 



The polynomials /, g, s are hereby defined by: 

f[rs\^l + l3irl +p2rs (25a) 
9[rs\^\rl + ^P2rs (25b) 

s[rs\^ ^r] + '^P2rs (25c) 

and the coefficients /3; and /14 are taken from Ref. |30] 
(where hi = ai~ hi(3^as in Table IV of Ref. |30]|: 

(3i = 1.158 1 3, 132 = 0.34455, K = -0.2942 (26) 

The dependence of the total capacitance C'^^^^ of a sym- 
metric capacitor on the electron density parameter r^, on 
the effective mass, and on the dielectric constant of the 
electrodes is presented in Figs. |4] and |5] for parameter 
sets with m*/m = 1 and 5, respectively. The contour 
lines of constant capacitance C^^^^ in Fig. [HJare qualita- 
tively different from those of the 2D capacitance (Fig. |3}, 
because the Thomas Fermi screening length with its func- 
tional dependence on m* and enters the expression for 
(7(30) _ Consequently, the transition to a negative capac- 
itance occurs at higher dielectric constants eos in the 3D 
electrodes than for the ideal 2D capacitor (cf. the scales 
of Eeff in Fig. [3] and Fig. |6]). For e^s = 1, however, the 
transition of a capacitor with 3D electrodes takes place 



at higher effective mass and larger carrier spacing (i.e. 
lower density) than for the 2D electrodes. 

We finally address the question, how these results com- 
pare to those of Ref. [THwhere e^s = 1 and = 1 have to 
be chosen.!^ Since we did not include the surface dipole 
layers, we cannot expect to find agreement concerning the 
effective capacitance length as defined by the inverse total 
capacitance. We therefore compare the dependencies of 
the kinetic term and of the combined exchange and cor- 
relation terms on Ts- We take the centre of mass position 
of the induced charge density from a recent publication 
(Table 1 in Ref. [59} in which the values with and without 
exchange and correlation are explicitly listed. The total 
values of the centre of mass position agree with those of 
Ref. [HI at rg = 2 and 4, except for a factor of 2 in Ref. [Til 
which has been suppressed in all subsequent publications: 

The inverse kinetic capacitance, which they obtain, in- 
creases similarly with as the one we are finding. The 
inverse capacitance associated with exchange and corre- 
lation energies is also negative in their data. The mag- 
nitude of the capacitance also increase with r^, however 
at a smaller rate. Our calculation based on a homoge- 
nous electron gas apparently overestimates the exchange 
effect. The overall trend of the contributions to the ca- 
pacitance reported in Ref. [59] agrees well with our find- 
ings, which seems surprising if the crude treatment of the 
surface inhomogeneity is considered.^^ 
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IV. CAPACITORS COMPRISING 
ELECTRODES WITH CORRELATED 
INHOMOGENEOUS ELECTRON SYSTEMS 

In a homogeneous electron gas, which resides in a con- 
tinuous medium, the compressibiHty of the interacting 
electrons is dominated by the kinetic and by the exchange 
contributions. In contrast, an electronic system that re- 
sides on a lattice may behave in a rather different man- 
ner, as epitomized by the Hubbard model in which the 
correlation term dominates the physics close to the Mott 
transition. Inter alia, it is the strong on-site interaction 
of charge carriers in the lattice models which on the one 
hand makes these models difficult to treat but on the 
other hand generates highly intriguing properties. Such 
strongly correlated systems may display a multitude of 
phase transitions, sometimes induced by small variations 
of the electronic density. But also, by altering tempera- 
ture, magnetic field, pressure, or epitaxial strain, phase 
transitions may be induced. Intermediate to strong elec- 
tronic correlations have been identified in systems such 
as the manganites, which display a colossal magnetoresis- 
tance effect, or such as the cuprates with their enigmatic 
pseudogap and superconducting states. 

Although for small variations of thermodynamic vari- 
ables the compressibility of the electronic system is often 
a dull thermodynamic quantity, close to phase transi- 
tions the compressibility can vary strongly. Correspond- 
ingly, the capacitance of an electrode built from such 
a material will reflect these variations in the electronic 
state — specifically the density dependence of the elec- 
tronic energy. It is pointed out that the compressibility 
is essentially proportional to the density of states of the 
interacting electronic system. Consequently, the capaci- 
tance depends on spectral weight transfer and formation 
of coherence peaks in the momentum integrated spec- 
tral function. Both of these phenomena are typical for 
strongly correlated electron systems and can, to a cer- 
tain extent, be controlled by temperature and magnetic 
field. It is therefore expected that they can produce ca- 
pacitive effects for electrodes. This prediction has to be 
investigated further. 

Whereas one does not expect the compressibility to 
vary appreciably close to the Mott transition in the single 
band Hubbard model, for the multi band case the situ- 
ation is markedly different: a large interorbital charge 
transfer with negative compressibility is feasible in a 
multiband model with at least one band being close to 
a Mott transition, as has been suggested by Liebsch for 
the insulator-metal transition in Sr-doped LaTiOg.ESln 
such multi-band (bulk) systems, the different orbitals act 
like the electrodes of a capacitor with negative compress- 
ibility at one "electrode" and zero distance between the 
"electrodes". Similar to these multiband correlated mate- 
rials, one may expect for mixed or intermediate valence 
systems a stronger dependence of the compressibilty on 
control parameters such as doping and temperature. 

A suppression of the density of states at the Fermi sur- 



face is known to occur in disordered correlated electronic 
systemd2II_ Al'tshuler and Aronov identified a square root 
suppression of the density of states at the Fermi level with 
decreasing temperature as a consequence of the inter- 
ference between inelastic and multiple, elastic scattering 
processes. Therefore we expect a corresponding decrease 
of the capacitance in case the surface of the electrodes 
consists of such a material. 



V. CAPACITORS WITH MULTIPLE STABLE 
OR METASTABLE STATES 

The energy functional of the capacitor is given by 
Eq. ([3}. The energy i? as a function of charge per 
area, Q/A, is identified from the dependence of E on 
the charge carrier density. For the model capacitor 
(Sec. Ill A I, for example, the energy per area is plotted as 
a function of charge density in Fig. |7] The curvature at 
Q/A ^ represents the inverse differential capacitance 
C^^{Q = 0) = dV/dQ which was calculated in Sec.|lII A 



Due to the negative capacitances of this model system, 
all electrons will accumulate on one of the two equivalent 
parallel plates, leave holes on the other plate, and thereby 



o 

o 

X 




Q/A xlO' [e/ae'] 



FIG. 7: Energy per area of a two-dimensional, parallel-plate 
capacitor plotted as a function of charge per area for several 
spacings of the plates. The energy has been calculated for 
the model capacitor described in Sec. [mX] The elect ronic 
systems of the electrodes form dilute and homogenous 2D 
electron gases. The charge carrier density (carriers per area) 
is W~'^/a% = 3.6 X lO^'^cm^^, the effective mass of the carri- 
ers is m* /m = 1, and the effective dielectric constant in the 
electrodes is ecff — 1 for all lines. The effective distance be- 
tween the plates d* /as = d/(e,. as) is 10.0, 8.0, 6.0, 5.0, and 
3.0 from top to bottom. 



11 




QIA xlO^ [elas^] 

FIG. 8: Sketch, illustrating several possible scenarios for 
the charge dependence of the energy of the two-dimensional, 
parallel-plate capacitor (compare Fig.[7|. The continuous red 
line corresponds to the red line of Fig. [Tjwith d* /as = 3.0. 
The dotted and dash-dotted lines are tentative extrapolations 
for high charge densities of the order of the carrier density 
10"^/a| = 3.6 X 10"cm"^. The latter lines have not been 
calculated, but rather present a sketch. The capacitor with 
four minima has different capacities in the metastable and 
stable states, as reflected by the different curvatures of the 
characteristics in the minima. Here, Ci and C2 refer to the 
inverse of the curvature at the minima which are marked by 
the vertical arrows. The capacitance of such a capacitor can 
therefore be switched by current pulses. 

charge the capacitoil^ (see the lower two curves in Fig. [tJ . 
The capacitor is characterized by multiple stable states, 
and may be switched between these by small bias pulses. 
For a capacitor with d* /as = d/ (e^ as) ~ 5.0 (the second 
curve from the buttom, in orange) and an area oi A = 
100 nm^, the energy difference to the zero charge state 
is approximately 0.5 eV, and the states will be stable 
aganst thermal activation. In this model, also capacitors 
with three stable or metastable states can be realized (the 
third curve from the buttom, in green), however with an 
energy difference of only about 50 meV. 

The model is not suited to precisely describe the 
multiple-state properties of real capacitors, because the 
nonlinear response of an inhomogeneous electron system 
will typically stabilize a charged state at charge densi- 
ties that are smaller than provided by the model. This 
scenario is illustrated in Fig. [s] (dotted curve). Phase 
transitions in electrodes with strongly correlated, inho- 
mogeneous electron systems might even produce addi- 
tional minima in the charge dependence of the energy 
functional (see Fig. [s] dash-dotted curve). The inves- 
tigation of such multistable charge states in capacitors 



that only use dielectric and not ferroelectric dielectrics is 
considered relevant for electronic devices such as memory 
devices or Qbits. 

Capacitors with negative differential capacitance gen- 
erate in the charged state electric fields in the dielectric 
and an electrical voltage between their electrodes. How- 
ever, the charged state is typically in thermal equilib- 
rium so that the electro-chemical potential vanishes and 
no charge flows through the leads by which the two elec- 
trodes are connected. It is noted that while such capac- 
itors are charged when being in one of the equilibrium 
states, these states are obviously stable states. The ca- 
pacitors do therefore not, as all other capacitors do, loose 
their charge with time by discharging. Treated as com- 
ponents in electronic circuits, capacitors with negative 
capacitances, such as the one for which the E{Q) char- 
acteristic is sketched by the continuous line of Fig. |8] 
will be characterized by the impedance Z = —i/{LuC) if 
operated close to the energy maximum. Thus, the fre- 
quency dependence of their impedance is the one of a 
standard capacitor. Because C is negative, however, the 
induced phase shift has the opposite sign of the phase 
shift of standard capacitors, so that the phase shift cor- 
responds to the phase shift of a standard inductance. In 




yxlO^ [Q/4neoaB] 

FIG. 9: Charge per area of a two-dimensional, parallel-plate 
capacitor plotted as a function of voltage between the plates. 
The charge per area Q{V) has been calculated for the model 
capacitor described in Sec. [Ill A| The electronic systems of 
the electrodes form dilute and homogenous 2D electron gases. 
The charge carrier density (carriers per area) is 10~^/a^ = 
3.6 X 10^"^cm~^, the effective mass of the carriers is m* /m = 
1, and the effective dielectric constant in the electrodes is 
eoff = 2.5. The effective distance between the plates d* /as = 
d/{eraB) is 10.0, where tr is the dielectric constant of the 
dielectric. The dotted line shows the Q(V^)-characteristic of 
the conventional, classical capacitor with Q = CgeomV'. 
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analogy, we suggest the possible realization of conduc- 
tors with negative inductances, the phase shift of which 
corresponds to the phase shift of standard capacitances. 

It is also amusing to consider the threshold behavior at 
a possible electrical breakdown of the dielectric. When 
approaching breakdown, localized charge carriers in the 
dielectric start to get released and move to the electrode 
of opposite charge. This charge redistribution shifts the 
electro-chemical potential so that an equal number of 
charge carriers flows through the leads. The original 
charge state is hereby stabilized and the release of charges 
does therefore not seem to change the electronic state 
in the electrodes. At the same time, the change of the 
charge distribution in the dielectric reduces the release 
of further carriers in the dielectric. The Coulomb-field of 
the altered charge distribution, however, affects the en- 
ergy of the electron systems of the plates and thereby in- 
directly alters the charge distribution on the electrodes. 
Should at breakdown a conducting channel be induced 
in the dielectric, the same argument applies, if the chan- 
nel is considered as an additional lead that connects the 
electrodes. The electronic state of the electrodes is not 
dependent on this channel and the electrodes will stay 
charged. 



VI. VOLTAGE DEPENDENCE OF THE 
CAPACITANCE 

In many applications it is the voltage and not the 
charge, with which the state of a capacitor is controlled. 
To calculate the C{V) characteristic we consider the Leg- 
endre transform £ of the energy E{Q) with respect to Q, 
and identify the total capacitance as 



CiV) 



Qjv) 

V 



(27) 
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FIG. 10: Capacitance C{V) = Q{V)/V of a two-dimensional, 
parallel-plate capacitor plotted as a function of voltage. The 
capacitance C{V) has been calculated for the model capacitor 
described in Sec. |III A[ The parameters are the same as in 
Fig. [9] The capacitance is normalized to Co, the capacitance 
of a classical capacitor of the same geometry. 




-2 2 

VxlO^ [Q/47I6oaB] 

FIG. 11: Capacitance and charge of the two-dimensional, 
parallel-plate capacitor of the model capacitor described in 
Sec. |III Aj plotted as a function of the voltage between the ca- 
pacitor plates. The charge carrier density (carriers per area) 
is /a% = 3.6x lO^'^cm"'^, the effective mass of the carriers 
is m* /m — 1, and the effective dielectric constant in the elec- 
trodes is Ecff ~ 1. The effective distance between the plates 
d* /as ~ d/{eraB) is 5.6. Top panel: voltage dependence 
of the capacitance C{V) = Q{V)/V. Lower panel: voltage 
dependence of the charge per area Q{V). In both panels, 
blue and red lines refer to the ranges of positive and neg- 
ative capacitances, respectively. The dotted line shows the 
characteristics of the conventional, classical capacitor with 



As example we take again the model system of Sec. |III A| 
with parallel, two-dimensional plates comprising the ca- 
pacitor. 

The calculations reveal that already in the regime of 
moderate effective distance of the capacitor plates, devia- 
tions from the conventional characteristics of a textbook 
capacitor are manifest: in Fig. |9]we display Q{V) which 
is expected to be linear with Q — CgcomV^ in the text- 
book case or for sufficiently large d* — djtr (green dot- 
ted line). For d* = lOa^ and = 2.5 the slope at zero 
voltage is steeper and Q(l^) deviates from a linear behav- 
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ior, most prominently for the strongly charged capacitor. 
The corresponding voltage dependent capacitance C{V) 
is shown in Fig. [lO] With increasing voltage the capac- 
itance increases because the capacitance is a function of 
the charge in the plates. The higher the charge in the 
plates the higher the capacitance which results in a no- 
ticeable upturn of C{V). 

This feedback of increasing capacitance with increasing 
charge may be so strong for nanoscopic distances that the 
conventional Q{V) and C{V) characteristic is completely 
modified (see Fig. [TT|. For example, for d* = S.Gob (and 
Ceff = 1) we operate in a regime where the capacitor 
self-charges once a sufficiently high charge has been ac- 
cumulated. In that case charge is transfered at vanishing 
energy expense (cf. Fig. [7) where V{Q) = dE/dQ = 
for a distinct, finite Q in the green, middle curve). The 



voltage is not zero in the fully charged state (Fig. Ill 



because this is not a thermodynamic state with a well- 
defined minimum in the energy E(Q) — as compared to 
the tentative capacitors of Fig. |8] which display minima 
in E{Q). 

The capacitance C{V) of the model capacitor (top 
panel of Fig. Ill may have multiple values, depending 



on the charge state. It is conceivable to switch between 
a positive value of C{V) (blue curve in Fig. Ill and a 



negative value (red curve in Fig. Ill, using charge pulses 



to switch the capacitor in and out of the "self-charged" 
state. 



VII. SPECIFIC EXAMPLES 

Magnetocapacitive effects have been investigated for 
several decades. A capacitive measurement, in which 
Landau levels in a 2D electron gas of a MOSFET were 
observed, was reported in Refs. I64|65I66I68|69I More re- 
cently, the density of states in the fractional quantum 
Hall regime has been investigated by magnetocapacitive 
measurements. A charge controlled transition of an 
electrode into a magnetic state would also yield large 
magnetocapacitive effects. Such a transition is conceiv- 
able in correlated electronic systems, for example, in the 
manganites. A positive magnetocapacitance has indeed 
been observed in Lag.ySro.sMnOa-^-titanate junction^Sl 
which contrasts to the negative magnetocapacitance of 
Pd-AlOaj-Al thin-film structures.''2lji Whereas the latter is 
the predicted behavior for a paramagnetic metal in which 
the Zeeman-split narrow d-bands cause the magnetoca- 
pacitance, the positive magnetoresistance has been re- 
lated to strong electronic correlations. 

Artificial materials that, for example, consist of lay- 
ered structures such as superlattices, of which the ca- 
pacitance is influenced by non-geometrical capacitances 
are predicted to be characterized by unusual properties. 
Of particular relevance are superlattices with a geometry 
chosen that th e ca pacitances between the layers are neg- 
ative (see Sec. Ill I so that the layers self-charge, and in 



netic phase transitions. Applied electric fields, on the 
one hand, change the carrier density of the magnetic lay- 
ers and therefore their magnetic properties. Application 
of magnetic fields, on the other hand, alters the correla- 
tion energies and therefore the capacitance and electric 
fields of the system. Thus, such artificial materials are 
candidates for strongly coupled, room temperature mul- 
tiferroics. 



A. Electrodes with 2D Electron Systems in 
Transition Metal Oxide Heterostructures 

The possibility to control the capacitance of capaci- 
tors by adding a thin, conducting sheet to an electrode 
opens the question how large a capacitance of a capacitor 
is expected that uses as an electrode a tw o-dim ensional 
electron gas generated at an oxide interfac^^^ED Indeed, 
capacitors based on two-dimensional electron gases in ox- 
ide heterostructures are expected to have unique and pos- 
sibly technically useful properties. 

The current model systems for such interfaces are 
systems like the n-doped LAO/STO or LVO/STO het- 
erostructures. The following discussion will refer to such 
interfaces, although future multilayers that do not imply 
sophisticated epitaxial growth may be preferable for ap- 
plications. For simplicity, we neglect strong correlation 
effects possibly present at such interfaces. 

As shown by Eqs. ( |8|9|12|13| , the capacitance (Eq. (|6j) 
of a capacitor that uses such an electron gas is a function 
of Ts, of the effective mass of the interface charge carriers 
m*, and of the effective dielectric constant ecs in the 
interface sheet. 

Up to now, such oxide interfaces have been prepared 
with sheet carrier densities of the order of ~ 5 x 10^^ — 
5 X 10^'*/cm^, corresponding to bare r^-values ranging 
from 50 to 5. As Fig. |3] shows for the example of a 
parallel-plate capacitor, that comprises a 1.5 nm thick 
Hf02 layer as dielectric (with — 23) and an electron 
gas with = 7.0 in each of the electrodes, the capaci- 
tances can be negative for realistic sets of to* and Ccff val- 
ues. Moreover, the capacitance is enhanced for a broad 
range of eofi and to*. 



which a part of the layers consist of materials with mag- 



B. Other Systems 

1/^kin may be strongly reduced at van Hove singular- 
ities. Conversely, the capacitance of condensators built 
with electrode materials such as graphene or other ma- 
terials with a small density of states, will be reduced be- 
cause then 1/Ckin dominates 1/Cgcom- In that case, the 
small density of states determines the charging of the ca- 
pacitance, and the charge is smaller than expected from 

C'goom • 

Recently, the quantum capacitance of coaxially gated 
carbon nanotube^^ has been calculated using a Green's 
function density functional tight-binding approach which 
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is based on a local density approximation (LDA) of the 
exchange correlation functional. This quantum capaci- 
tance includes the terms Ckin and Cx which are evaluated 
for the quasi one-dimensional nanotube. The behavior 
of this capacitance (see Fig. 3 in Ref. [71} is consistent 
with our results for the capacitance of 2D or 3D capaci- 
tors. Specifically, in the regime of small carrier densities 
overscreening of the gate field was observed which corre- 
sponds to a negative capacitance of the carbon nanotube. 
The capacitance of carbon nanotube has been measured 
by a group from Cornell UniversitjEH and, in fact, the 
measurements suggest the existence of a negative contri- 
bution to the total capacitance. 

A negative electronic compressibility is also expected 
close to structural transitions. Commonly, this signifies 
that the electronic system is in a phase separated state 
or undergoes charge disproportionation. In many cases, 
these situations are associated with a metal-insulator 
transition. If structural transitions can achieve a func- 
tional quality for a controlled electronic transition in the 
electrodes of a capacitor has to be investigated. 



material of the electrodes. It suffices to introduce thin 
conducting sheets of a compound (see Fig. 
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3 ) with the 
ter the ca- 



appropriate electron energies to drastically a 
pacitance. 

Therefore, besides optimizing the thickness and the di- 
electric properties of the dielectric layer, choosing the 
appropriate materials of the electrode surfaces provides 
an independent access to optimizing the capacitance of 
capacitors. Materials of interest include materials with 
small electron densities, such as interfaces in oxide het- 
erostructures, low-electron density metals, materials with 
strong electronic correlations such as transition metal ox- 
ides, or graphene-based systems. 

Because the capacitance of such capacitors is controlled 
by the energy of the electron system of the plates rather 
than by the field energy, the standard parallel-plate con- 
figuration is not needed to obtain large capacitances or 
field strengths (for example, in Fig. 12 1 all blue surfaces 
contribute to the capacitance 1/C — 1/Cgoom in Eq. (|6|). 
Rather, the surface areas of the electrodes are important. 
We note that the desired modification of the capacitance 
does not require that the electrodes consist entirely of 



VIII. SUMMARY, CONCLUSIONS AND 
PERSPECTIVES 



Considering the complete energy of capacitors, which 
includes the Coulomb energy and the energy of the elec- 
tron systems of the capacitor plates, we have evaluated 
the capacitance of capacitors. We find that in general 
the capacitance is given by Eq. ^ and not by Eq. ([TJ. 

As shown, the capacitance of a capacitor is described 
by the series connection of the Coulomb capacitance, the 
kinetic capacitance of the electrons, the exchange capac- 
itance, the correlation capacitance, and the capacitance 
which results from the interaction of electrons with non 
charge-carrier degrees of freedom. Because for a free elec- 
tron gas the exchange capacitance and the correlation 
capacitance are negative and the electron-phonon inter- 
action may also result in a negative value, the total ca- 
pacitance may be smaller or larger than the Coulomb 
capacitance, and may even reach negative values. 

Our analysis shows that the capacitance is dominated 
by the field capacitance for the case that the field energy 
exceeds the energies of the electron system, but that the 
electron energies are significant if the thickness of the 
dielectric layer divided by its dielectric constant is com- 
parable to the Bohr radius. This applies independent of 
the lateral size of the capacitor. The gate capacitances 
of modern MOSFETs are reaching this regime. 

Eq. (|6j is of general importance to calculate the capac- 
itances of capacitors. Due to their broad applicability, we 
have specifically provided the detailed calculations of the 
capacitance of parallel-plate capacitors with electrodes 
that contain two-dimensional and three-dimensional, ho- 
mogenous electron systems. 

It is concluded that the capacitance of capacitors can 
be significantly enhanced or decreased by altering the 
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plate 1 
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plate 2 



FIG. 12: Sketch, illustrating options to integrate materials 
into capacitors to achieve capacitances that deviate strongly 
from the geometric capacitance. Panel (a) sketches a situ- 
ation where electrodes may cause a strong variaton of the 
capacitance due to their contributions in Cj^^^, C^^, C^^ or 



c„- 



so that a large total capacitance is obtained with- 



out using a parallel-plate configuration.. Panel (b) sketches 
a thin, conducting sheet (booster layer) of a compound with 
appropriate electron energies integrated into a standard ca- 
pacitor. This sheet can considerably enhance or reduce the 
capacitance. The tickness of the sheet may be as small as the 
electric screening length in the compound. 
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FIG. 13: Sketch to illustrate the possibility of electric field 
enhancement due to the integration of a material with high 
capacitance into a plate of a capacitor. The material is in- 
troduced in a partial coating (blue stripes) of the capacitor 
plates (dark gray rectangles). 



strongly correlated compounds. It is sufficient to coat 
the active surface of an electrode with a thin layer of a 
correlated material, because the net charges are repelled 
from the interiors of the electrodes. 

We emphasize that eofi is a property of the electrode 
materials and, specifically, it relates to the properties of 
the materials at their surfaces, ecs therefore can be tuned 
by coating the electrodes with thin layers of desired ma- 
terials and then is given by the dielectric constant of the 
coating material. It accounts for the polarizability of the 
ionic cores and also includes the short-range screening 
processes from excitations in the adjacent material lay- 
ers, but not the contributions from the electron gas itself. 
The effective dielectric constant ees is, beyond linear re- 
sponse, also dependent on the charge density. Such a 
dependence induces feedback effects which alter the ca- 
pacitance. While one may consider to estimate from 
an LDA evaluation for the specific heterostructure, there 
is no simple scheme to calculate Coff. Lacking such a 
simple scheme, we calculated the capacitance of a simple 
"model capacitor" for several values of eog. 

Using optimized capacitor configurations, it is in prin- 
ciple possible to build capacitors that are characterized 
by a negative capacitance for zero applied voltage. Such 
capacitors charge themselves to generate finite electric 
fields and, although they do not contain any ferroelectric 
component, show parallels to ferroelectric devices. It is 
obvious, that for device applications novel possibilities 
are opened by adding truly ferroelectric layers to struc- 
tures with unconventional capacities. Intriguing charac- 
teristics arise if the capacitances are negative. 

In the case of negative compressibility, adding or re- 
moving an electron from an electrode involves less energy 
than the standard work function or electron affinity. It 
is intriguing to regard the behavior of the electrons in a 
single electrode or on the two electrodes of a capacitor, in 



case the compressibility of at least one electron system 
is negative. In a capacitor that self charges, for exam- 
ple, into a positively and negatively charged plate, the 
electrons and the holes in the two plates interact across 
the dielectric by their Coulomb field. It seems possible 
that electron-hole pairs result that form a coherent state, 
and, for example modify pairing in a superconducting 
electrode. Such effects may be particularly useful in case 
the plate capacitor is designed as a transmission line. We 
propose that such a superconducting state may even ex- 
ist if none of the electrodes is superconducting by itself. 
Even for electron systems with positive compressibilities, 
the Coulomb energies caused by the presence of a second 
plate will lower or enhance the Tc of a superconducting 
electrode. Therefore, modifying Coulomb energies by us- 
ing capacitor-type configurations, the electrodes of which 
may even consist of the same materials, offers new pos- 
sibilities to realize superconductors. 

Also, because the capacitance can be controlled by the 
electronic system of the plates, which may be strongly 
correlated electron systems that are sensitive to mag- 
netic fields, materials may be designed from capacitor- 
derived multilayer configurations that yield large multi- 
ferroic coupling at room temperature. 

Capacitors for which the non-geometrical capacitances 
are important, and in which one plate has a non- 
homogeneous composition, are characterized by inhomo- 
geneous and unusual electric field distributions. The 
charge distribution of such an electrode is modified and 
differs from the charge distribution generated in the cor- 
responding standard configuration. Consequently, in the 
presence of an electric field, the distribution of the field 
lines is accordingly altered. As illustrated in Fig. 13 
electric field lines can, for example, be focused onto de- 
sired areas to extract charges. Or, where desirable, areas 
can be shielded from electric fields, which may be impor- 
tant to prevent electric breakthrough phenomena. Thus, 
inhomogeneous electrodes allow to generate desirable dis- 
tributions of charge densities and electric fields without 
altering the overall geometrical configuration of a device. 
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